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Abstract. We study ways to expedite Yates's algorithm for computing the zeta and 
Moebius transforms of a function defined on the subset lattice. We develop a trimmed 
variant of Moebius inversion that proceeds point by point, finishing the calculation at 
a subset before considering its supersets. For an n-element universe U and a family 5" 
of its subsets, trimmed Moebius inversion allows us to compute the number of packings, 
coverings, and partitions of U with k sets from J in time within a polynomial factor (in 
n) of the number of supersets of the members of 3. 

Relying on an intersection theorem of Chung et al. (1986) to bound the sizes of set 
families, we apply these ideas to well-studied combinatorial optimisation problems on 
graphs of maximum degree A. In particular, we show how to compute the Domatic Number 
in time within a polynomial factor of (2 A+1 — 2) n,/ ' A+1 ' and the Chromatic Number in 
time within a polynomial factor of (2 A+1 — A — 1)™/( A + 1 ). For any constant A, these 
bounds are 0((2 — e) n ) for e > independent of the number of vertices n. 



1. Introduction 

Yates's algorithm from 1937 is a kind of fast Fourier transform that computes for a 
function / : {0, l} n — » R and another function v : {0, 1} x {0, 1} — » R the values 

f(xi, . . . , x n ) = ^2 v(x 1 ,yi)---v(x n ,y n )f(y 1 ,...,y n ). (1.1) 

m,-,VnE{o,i} 

simultaneously for all X = (x\, . . . , x n ) € {0, l} n using only 0(2 n n) operations, instead of 
the obvious 0(4 n n). The algorithm is textbook material in many sciences. Yet, though it 
appears in Knuth |X3|, §3.2], it has received little attention in combinatorial optimisation. 

Recently, the authors [21 [3] used Yates's algorithm in combination with Moebius inver- 
sion to give algorithms for a number of canonical combinatorial optimisation problems such 
as Chromatic Number and Domatic Number in n-vertex graphs, and n-terminal Minimum 
Steiner Tree, in running times within a polynomial factor of 2 n . 
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From the way it is normally stated, Yates's algorithm seems to face an inherent 2 n 
lower bound, up to a polynomial factor, and it also seems to be oblivious to the structural 
properties of the transform it computes. 

The motivation of the present investigation is to expedite the running time of Yates's 
algorithm for certain structures so as to get running times with a dominating factor of the 
form (2 — e) n . From the perspective of running times alone, our improvements are modest at 
best, but apart from providing evidence that the aesthetically appealing 2 n bound from [3] 
can be beaten, the combinatorial framework we present seems to be new and may present 
a fruitful direction for exact exponential time algorithms. 

1.1. Results 

In a graph G = (V, E), a set D C V of vertices is dominating if every vertex not in D 
has at least one neighbour in D. The domatic number of G is the largest k for which V can 
be partitioned in to k dominating sets. We show how to compute the domatic number of 
an n-vertex graph with maximum degree A in time 

0*((2 A+1 -2) n /( A+1 )) ; 

the O* notation suppresses factors that are polynomial in n. For constant A, this bound is 
always better than 2 ra , though not by much: 



A 


3 


4 


5 


6 


7 


8 


(2 A+1 - 2) 1/(A+1) 


1.9344 


1.9744 


1.9895 


1.9956 


1.9981 


1.9992 • • • 



The chromatic number of a graph is the minimum k for which the vertex set can be 
covered with k independent sets; a set / C V is independent if no two vertices in / are 
neighbours. We show how to compute the chromatic number of an n-vertex graph with 
maximum degree A in time 

0*((2 A + 1 _A-l) n /( A+1 )). 
This is slightly faster than for Domatic Number: 



A 


3 4 5 


6 


7 8 


( 2 A+1 - A - 1)V(A+1) 


1.8613 1.9332 1.9675 


1.9840 


1.9921 1.9961 ••■ 



One notes that even for moderate A, the improvement over 2 n is minute. Moreover, 
the colouring results for A < 5 are not even the best known: by Brooks's Theorem [5], 
the chromatic number of a connected graph is bounded by its maximum degree unless 
the graph is complete or an odd cycle, both of which are easily recognised. It remains to 
decide if the chromatic number is 3, 4, or 5, and with algorithms from the literature, 3- 
and 4-colourability can be decided in time 0(1.3289 n ) pQ and 0(1. 7504") [6], respectively. 
However, this approach does stop at A = 5, since we know no o(2 n ) algorithm for 5- 
colourability. Other approaches for colouring low-degree graphs are known via pathwidth: 
given a path decomposition of width w the /c-colourability can be decided in time k w n°^ 
|llj : for 6-regular graphs one can find a decomposition with w < n(23 + e)/45 for any e > 
and sufficiently large n [TT], and for graphs with m edges one can find w < m/5.769 + 
O(logn) [12J. However, even these pathwidth based bounds fall short when k > 5 — we are 
not aware of any previous o(2 n ) algorithm. 

For the general case, it took 30 years and many papers to improve the constant in 
the bound for Chromatic Number from 2.4423 [14] via 2.4151 [9], 2.4023 [6], 2.3236 [2], 
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Figure 1: Trimmed evaluation. Originally, Yates's algorithm considers 
the entire subset lattice (left). We trim the evalation from below by 
considering only the supersets of 'interesting' points (middle), and from 
above by abandoning computation when we reach certain points (right). 

to 2 [4|, and a similar (if less glacial) story can be told for the Domatic Number. None 
of these approaches was sensitive to the density of the graph. Moreover, what interests us 
here is not so much the size of the constant, but the fact that it is less than 2, dispelling 
the tempting hypothesis that 2 n should be a 'difficult to beat' bound for computing the 
Chromatic Number for sparse graphs. In £J5] we present some tailor-made variants for 
which the running time improvement from applying the ideas of the present paper are more 
striking. 

Chromatic Number and Domatic Number are special cases of set partition problems, 
where the objective is to partition an n-element set U (here, the vertices of a graph) with 
members of a given family 3~ of its subsets (here, the independent or dominating sets of the 
graph). In full generality, we show how to compute the covering, packing, and partition 
numbers of (U, H) in time within a polynomial factor of 

\{T CU : there exists an S G J such that S CT}\, (1.2) 

the number of supersets of the members of 3 r . In the worst case, this bound is not better 
than 2 n , and the combinatorial challenge in applying the present ideas is to find good bounds 
on the above expression. 

1.2. Techniques 

The main technical contribution in this paper, sketched in Figure 1, is that Yates's 
algorithm can, for certain natural choices of v : {0, 1} x {0, 1} — > R, be trimmed by 
considering in a bottom-up fashion only those X G {0, l} n that we are actually interested 
in, for example those X for which f(X) ^ and their supersets. (We will understand 
X as a subset of {1, . . . ,n} whenever this is convenient.) Among the transforms that are 
amenable to trimming are the zeta and Moebius transforms on the subset lattice. 

We use the trimmed algorithms for zeta and Moebius transforms to expedite Moebius 
inversion, a generalisation of the principle of inclusion-exclusion, which allows us to compute 
the cover, packing, and partition numbers. The fact that these numbers can be computed 
via Moebius inversion was already used in [21 S] , and those parts of the present paper 
contain little that is new, except for a somewhat more explicit and streamlined presentation 
in the framework of partial order theory. 

The fact that we can evaluate both the zeta and Moebius transforms pointwise in such 
a way that we are done with X before we proceed to Y for every Y D X also enables us to 
further trim computations from what is outlined above. For instance, if we seek a minimum 
set partition of sets from a family 3~ of subsets of U, then it suffices to find the minimum 
partition of all X such that U \ X = S for some S G IT. In particular, we need not consider 
how many sets it takes to partition X for X's large enough for U \ X not to contain any 
set from 3. 

The main combinatorial contribution in this paper is that if 5F is the family of maximal 
independent sets, or the family of dominating sets in a graph, then we show how to bound 
(]1.2p in terms of the maximum degree A using an intersection theorem of Chung et al. 
[8] that goes back to Shearer's Entropy Lemma. For this we merely need to observe that 
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the intersection of 3 and the closed neighbourhoods of the input graph excludes certain 
configurations. 

In summary, via (|1.2p the task of bounding the running time for (say) Domatic Number 
reduces to a combinatorial statement about the intersections of certain families of sets. 

Notation. Yates's algorithm operates on the lattice of subsets of an n-element universe U, 
and we find it convenient to work with notation established in partial order theory. 

For a family 3 of subsets of U, let min^F (respectively, max 3") denote the family of 
minimal (respectively, maximal) elements of 3 with respect to subset inclusion. The upper 
closure (sometimes called up-set or filter) of 3 is defined as 

13 = { T C U : there exists an S E 3 such that S C T } . 

For a function / defined on subsets of U, the support of / is defined as 

supp(/) = {AC[/:/(A)^0}. 

For a graph G, we let CD denote the family of dominating sets of G and J the family 
of independent sets of G. Also, for a subset W C V of vertices, we let G[W] denote the 
subgraph induced by W. For a proposition P, we use Iverson's bracket notation [P] to 
mean 1 if P is true and otherwise. 

2. Trimmed Moebius Inversion 

For a family 3 of sets from {0, l} n and a set A E {0, l} n we will consider /c-tuples 
(Si, . . . , Sk) with Si £ 3 and Sj C X. Such a tuple is disjoint if Sj x n Sj 2 =0 for all 
l<ii<i2<&> and covering if Si U ■ • ■ U Sfc = X. From these concepts we define for fixed 
k 

(1) the cover number c(X), viz. the number of covering tuples, 

(2) the packing number p(X), viz. the number of disjoint tuples, 

(3) the partition number or disjoint cover number d(X), viz. the number of tuples that 
are both disjoint and covering. 

In this section we show how to compute these numbers in time H^n ^, rather than 
2" n O(i) ag j n [3j gj. The algorithms are concise but somewhat involved, and we choose 
to present them here starting with an explanation of Yates's algorithm. Thus, the first 
two subsections are primarily expository and aim to establish the new ingredients in our 
algorithms. 

At the heart of our algorithms lie two transforms of functions / : {0, l} n — > R on the 
subset lattice. The zeta transform f( is defined for all X E {0, l} n by 

(/0(A) = £ /(F) . (2.1) 

YCX 

(The notation /£ can be read either as a formal operator or as a product of the 2 n - 
dimensional vector / and the matrix Q with entries Qyx = [F C A].) The Moebius transform 
ffi is defined for all A E {0, l} n by 

{f MX) = £(-l)l*\ y l/(F). (2.2) 

YCX 
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These transforms are each other's inverse in the sense that / = = f[i(, a fundamental 
combinatorial principle called Moebius inversion. We can (just barely) draw an example 
in four dimensions for a function / given by /({4}) = /({1,2, 4}) = 1, /({1,3}) = 2 and 
f(X) = otherwise: 




Another example that we will use later is the connection between the packing number and 
the disjoint cover number, 

p = d(, (2.3) 

which is easy to verify: By definition, 

(O(x) = £d(y). 

YCX 

Every disjoint fc-tuple (Si, . . . , Sk) with Si U • • • U Sk Q X appears once on the right hand 
side, namely for Y = S% U • • • U Sk, so this expression equals the packing number p(X). 

2.1. Yates's algorithm 

Yates's algorithm [T7] expects the transform in the form of a function v : {0, 1} x 
{0, 1} — > R and computes the transformed values 

f(X) = Yl v{xi,yi)---v{x n ,y n )f{Y). (2.4) 
Ye{o,i} n 

simultaneously for all X € {0, l} n . Here, we let {x\, . . . ,x n ) and (yi, . . . ,y n ) denote the 
binary representations (or, 'incidence vectors') of X and Y, so Xj = [j € X] and yj = [j € 
Y]. To obtain (|2.1|) set v(x,y) = [y < x] and to obtain (|2.2|) set v(x,y) = [y < x](—l) x ~ y . 

The direct evaluation of (|2.4f) would take 2 n evaluations of / for each X, for a total of 
0{2 n 2 n n) = 0(4 n n) operations. The zeta and Moebius transforms depend only on Y C X, 
so they would require only Ylx 2'^' = X^o<i<n (?)^' = ^™ evaluations. Yates's algorithm is 
faster still and computes the general form in 0{2 n n) operations: 

Algorithm Y. {Yates's algorithm.) Computes f(X) defined in §XZ$ for all X G {0, 1}™ given f(Y) 
for all Y e {0, 1}™ and v(x,y) for all x,y G {0, 1}. 

Yl: For each X G {0, 1}", set g (X) =f(X). 
Y2: For each j = 1, . . . , n and X G {0, 1}", set 

9j (X) = v(\j G X],0) 9j ^(X \ {j}) + v(\j G X], l)9i-i{X U {j}) . 

Y3: Output g n . 

The intuition is to compute f(X) 'coordinate-wise' by fixing fewer and fewer bits of X 
in the sense that, for j = 1, . . . , n, 

9i( X ) = ^ v(x 1 ,y 1 )---v(x j ,y j )f(y 1 ,...,y j ,x j+1 ,...,x n ). (2.5) 

yi,...,j/ J e{0,l} 
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Indeed, the correctness proof is a straightforward verification (by induction) of the above 
expression. 

2.2. Trimmed pointwise evaluation 

To set the stage for our present contributions, observe that both the zeta and Moebius 
transforms 'grow upwards' in the subset lattice in the sense that supp(/C), supp(//x) C 
|supp(/). Thus, in evaluating the two transforms, one ought to be able to trim off redundant 
parts of the lattice and work only with lattice points in |supp(/). 

We would naturally like trimmed evaluation to occur in 0(||supp(/)|n) operations, in 
the spirit of Algorithm Y. However, to obtain the values at X in Step Y2 of Algorithm Y, 
at first sight it appears that we must both 'look up' (at X U {j}) and 'look down' (at 
X\{j}). Fortunately, it suffices to only 'look down'. Indeed, for the zeta transform, setting 
v(x,y) = [y < x] and simplifying Step Y2 yields 

9j (X) = [j G X] 9j ^(X \ {j}) + (X) . (2.6) 

For the Moebius transform, setting v(x,y) = [y < x](—l) x ~ v and simplifying yields 

9j (X) = -[j £ X]gj-i(X \ {j}) +g,- 1 {X) . (2.7) 

Furthermore, it is not necessary to look 'too far' down: for both transforms it is immediate 
from ([23D that 

gj(X) = holds for all X <£ |supp(/) and j = 0, . . . , n . (2.8) 

In what follows we tacitly employ (|2.8p to limit the scope of (I2.6P and (12. 7|) to |supp(/). 

The next observation is that the lattice points in T SU PP(/) can be evaluated in order of 
their rank, using sets £(?") containing the points of rank r. Initially, the sets £(r) contain 
only supp(/), but we add elements from T SU PP(/) as we g° along. These observations result 
in the following algorithm for evaluating the zeta transform; the algorithm for evaluating 
the Moebius transform is obtained by replacing (|2.6I) in Step Z3 with (12. 7ft . 

Algorithm Z. ( Trimmed pointwise fast zeta transform.) Computes the nonzero part of /C given the 
nonzero part of/. The algorithm maintains n+1 families £(0), ■ ■ • ,&(n) of subsets X e {0, 1}™; £(r) 
contains only sets of size r. We compute auxiliary values gj(X) for all 1 < j < n and X G |supp(/); 
it holds that g n (X) = (fQ{X). 

Zl: For each X £ supp(/), insert X into £(|X|). Set the current rank r = 0. 
Z2: Select any X E £j(r) and remove it from £(r). 
Z3: Set g (X) = f{X). For each j = 1, . . . , n, set 

9j (X) = [j e X] gj -x(X \ {j}) + g 3 ^{X) . 

[At this point g n (X) = (f()(X).} 
Z4: If g n {X) ^ 0, then output X and g n {X). 
Z5: For each j £ X, insert X U {j} into L(r + 1). 

Z6: If iL(r) is empty then increment r < n until £(r) is nonempty; terminate if r = n and L(n) 

is empty. 
Z7: Go to Z2. 

Observe that the evaluation at X is complete once Step Z3 terminates, which enables 
further trimming of the lattice 'from above' in case the values at lattice points with higher 
rank are not required. 
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By symmetry, the present ideas work just as well for transforms that 'grow downwards', 
in which case one needs to 'look up'. However, they do not work for transforms that grow 
in both directions, such as the Walsh-Hadamard transform. 

In the applications that now follow, / will always be the indicator function of a family 
3". In this case having supp(/) quickly available translates to 3~ being efficiently listable; 
for example, with polynomial delay. 

2.3. Covers 

The easiest application of the trimmed Moebius inversion computes for each X E ]3 
the cover number c(X). This is a particularly straightforward function of the zeta transform 
of the indicator function /: simply raise each element of to the A;th power and transform 
the result back using \i. To see this, observe that both sides of the equation 

(c()(Y) = ((fC)(Y)) k (2.9) 

count the number of ways to choose A:-tuples (Si, . . . ,S k ) with Si C Y and Si G 3\ By 
Moebius inversion, we can recover c by applying fi to both sides of f|2.9[) . 

Algorithm C. (Cover number.) Computes c(X) for all X G j£F given £F. The sets £(r) and auxiliary 
values gj{X) are as in Algorithm Z; also required are auxiliary values hj(X) for Moebius transform. 

CI: For each leJ, insert X into £(|X|). Set the current rank r = 0. 

C2: Select any X E £(r) and remove it from £(r). 

C3: [Zeta transform.] Set go(X) = [IeJ]. For each j = 1, . . . ,n, set 

9j (X) = [j G X]g^ 1 (X \ {j}) + g 3 ^{X) . 

[At this point it holds that g n (X) = (f()(X).] 
C4: [Evaluate zeta transform of c(X).] Set ho(X) = g n (X) k . 
C5: [Moebius transform.] For each j — 1, . . . ,n, set 

h 3 (x) = -\j e x]h 3 ^{x \ {j}) + h^x) . 

C6: Output X and h n (X). 

C7: For each j X, insert X U {j} into H(r + 1). 

C8: If £j(r) is empty, then increment r < n until &(r) is nonempty; terminate if r = n and L(n) 

is empty. 
C9: Go to C2. 

2.4. Partitions 

What makes the partition problem slightly less transparent is the fact that we need to 
use dynamic programming to assemble partitions from sets with different ranks. To this 
end, we need to compute for each rank s the 'ranked zeta transform' 

(/C (s) )P0 = £ /CO- 

YCX,\Y\=s 

For rank s, consider the number S S \Y) of tuples (Si, . . . , S^) with Si E 3~, Sj C Y, 
Si U • • • U S k = Y and [Si| + • • • + \S k \ = s. Then d(Y) = d^ Y \\Y). Furthermore, the 
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zeta-transform (d^ s \)(X) counts the number of ways to choose (Si, . . . , Sk) with Si C X, 
Si G y, and \Si\ + • • • + \Sk\ = s. Another way to count the exact same quantity is 

k 

q(k,s,X) = Yl n^X*)- (2-10) 

siH hs k =s i=l 

Thus we can recover d&(Y) from q(k,s,X) by Moebius inversion. 

As it stands, (|2.10j) is time-consuming to evaluate even given all the ranked zeta trans- 
forms, but we can compute it efficiently using dynamic programming based on the recurrence 



ELo q(k-i,s-t, x)(f(W)(x), if k > l , 

(fC (s) )(X), ]£k = l. 



q(k,s,X) -- 
This happens in Step D4. 

Algorithm D. (Disjoint cover number.) Computes d(X) for all X € given 3\ The sets H(r) are 
as in Algorithm Z; we also need auxiliary values gj(X) and hi (X) for all X 6 "f J", 1 < j < n, and 
< s < n; it holds that g { n s) {X) = (f(^)(X) and h ( n\x) = d^(X). 

Dl: For each X £ 7, insert X into Set the current rank r = 0. 

D2: Select any X € £(»") and remove it from L(r). 

D3: [Ranked zeta transform.] For each s = 0, . . . , n, set c^ s) (X) = [Ie = s]. For each 

j = 1, . . . , n and s = 0, . . . , n, set 

5 j s) (X) = [.? e Xlffj'A^ \ {j}) + g^(X) . 



[At this point it holds that g { n s) (X) = (/C (s) )(A) for all < s < 



n . 



D4: [Evaluate zeta transform of dS s \] For each s = 0, ...,n, set q(l, s) = gn(X). For each 
i = 2, . . . , k and s = 0, . . . , n, set q(i, s) — J2t=o ?(* — 1j s — (^0- 

D5: [Ranked Moebius transform.] For each s = 0, . .. ,n, set Hq S \x) = q(k, s). For each j = 
1, . . . , n and s = 0, . . . , n, set 

h<f\x) = -\j e X]h^(X \ {j}) + h^(X) . 

[At this point it holds that h$(X) = d (s \X) for all < s < n] 
D6: Output X and h ( ,l Xl} (X). 
D7: For each j £ X, insert X U {j} into £(r + 1). 

D8: If £(r) is empty, then increment r < n until fL(r) is nonempty; terminate if r = n and £(n) 

is empty. 
D9: Go to D2. 



2.5. Packings 



According to (|2.3p . to compute p(X) it suffices to zeta-transform the partition number. 
This amounts to running Algorithm Z after Algorithm D. (For a different approach, see [3].) 
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3. Applications 

3.1. The number of dominating sets in sparse graphs 

This section is purely combinatorial. Let D denote the dominating sets of a graph. A 
complete graph has 2 n — 1 dominating sets, and sparse graphs can have almost as many: 
the n-star graph has 2 n_1 dominating sets and average degree less than 2. Thus we ask how 
large |D| can be for graphs with bounded maximum degree. An easy example is provided 
by the disjoint union of complete graphs of order A + 1: every vertex subset that includes 
at least one vertex from each component is dominating, so |D| = (2 A+1 — l) n /( A+1 ). We 
shall show that this is in fact the largest possible D for graphs of maximum degree A. Our 
analysis is based on the following intersection theorem. 

Lemma 3.1 (Chung et al. [8j). Let U be a finite set with subsets P\, . . . , P m such that 
every u G U is contained in at least 5 subsets. Let J be a family of subsets of U . For each 
1 < £ < m, define the projections 3"f = {FnPf : F 6 3"}. Then 

m 

m*<ni3ii- 

1=1 

Theorem 3.2. The number of dominating sets of an n-vertex graph with maximum degree 
A is at most (2 A+1 - 1)»/( A + 1 ). 

Proof. Let G = (V, E) be a graph with \V\ = n and maximum degree A. For each v G V, 
let A v be the closed neighbourhood around vertex v, 

A v = {v}U{u G V : uv G E} . (3.1) 

Next, for each u G V with degree d(u) < A, add u to A — d{u) of the sets A v not already 
containing u (it does not matter which). Let a v = \A V \ and note that a v = (A + l)n. 

We want to apply Lemma 13.11 To this end, let U = V and m = n. By construction, 
every u G V belongs to exactly S = A + 1 subsets A v . To get a nontrivial bound on T> 
we need to bound the size of D v = { D n A v : D G D }. Every D n A v is one of the 2 av 
subsets of A v , but none of the D n A v can be the empty set, because either v or one of its 
neighbours must belong to the dominating set D. Thus \D V \ < 2 av — 1. By Lemma 13-H we 
have 

\T)\ A+1 <]J{2 a - -1) . (3.2) 

V 

Since x i— > log (2 X — 1) is concave, Jensen's inequality gives 

- V log (2 a " - 1) < log (2^ a "/ n - 1) = log (2 A+1 - 1) . 
n ^ 

V 

Taking exponentials and combining with (|3,2p gives |D| A+1 < (2 A+1 — l) n . ■ 

3.2. Domatic Number 

We first observe that a graph can be packed with k dominating sets if and only if it 
can be packed with k minimal dominating sets, so we can consider /c-packings from min T> 
instead of T). This has the advantage that minD can be listed faster than 2 n . 

Lemma 3.3 (Fomin et al. [10J). Any n-vertex graph has at most O*(1.7170 n ) minimal 
dominating sets, and they can be listed within that time bound. 
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Theorem 3.4. For an n-vertex graph G with maximum degree A we can decide in time 



whether G admits a packing with k dominating sets. 

Proof. We use Algorithm D with 3~ = min T>. By the above lemma, we can complete Step Dl 
in time O*(1.7170 n ). The rest of the algorithm requires time 0*(| jminD|). Since every 
superset of a dominating set is itself dominating, jminD is a sub-family of T> (in fact, it is 
exactly T>), so Theorem 13.21 bounds the total running time by 



We can do slightly better if we modify Algorithm D in Step D7 to insert X U {j} only 
if it excludes at least one vertex for each closed neighbourhood. Put otherwise, we insert 
X U {j} only if the set V \ (X U {j}) dominates the graph G. The graph then has Domatic 
Number at least k+ 1 if and only if the algorithm reports some X for which d{X) is nonzero. 
The running time can again be bounded as in Theorem 13.21 but now D n A v can neither be 
the empty set, nor be equal to A v . Thus the application of Lemma 13. II can be strengthened 
to yield the claimed result. ■ 

3.3. Chromatic Number 

Our first argument for Chromatic Number is similar; we give a stronger and slightly 
more complicated argument in §3.41 

We consider the independent sets J of a graph. An independent set is not necessarily 
dominating, but it is easy to see that a maximal independent set is dominating. Moreover, 
the Moon-Moser bound tells us they are few, and Tsukiyama et al. tell us how to list them 
with polynomial delay: 

Lemma 3.5 (Moon and Moser [15J; Tsukiyama et al. [16J). Any n-vertex graph has at 
most O* (1.4423™) maximal independent sets, and they can be listed within that bound. 

Theorem 3.6. For an n-vertex graph G with maximum degree A we can decide in time 



whether G admits a covering with k independent sets. 

Proof. It is easy to see that G can be covered with k independent sets if and only if it can 
be covered with k maximal independent sets, so we will use Algorithm C on max J. Step CI 
is completed in time 0*(1.4423 n ), and the rest of the algorithm considers only the points 
in jmaxJ, which all belong to T). Again, Theorem 13.21 bounds the total running time. ■ 

3.4. Chromatic Number via bipartite subgraphs 

We can do somewhat better by considering the family S of vertex sets of induced 
bipartite subgraphs, that is, the family of sets B C V for which the induced subgraph G[B] 
is bipartite. As before, the literature provides us with a nontrivial listing algorithm: 

Lemma 3.7 (Byskov and Eppstein [7]). Any n-vertex graph has at most 0*(1.7724 n ) 
maximal induced bipartite subgraphs, and they can be listed within that bound. 

The family max 2 is more than just dominating, which allows us to use Lemma 13. II in 
a stronger way. 





0*((2 A+1 



l)n/(A+l)\ 



TRIMMED MOEBIUS INVERSION AND GRAPHS OF BOUNDED DEGREE 
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Theorem 3.8. For an n-vertex graph of maximum degree A it holds that 

| "[max £ | < (2 A+1 - A - 1)™/( A+1 ) . 

Proof. Let G = (V,E) be a graph with \V\ = n and maximum degree A. Let 3~ = f max 23. 
Let A v be as in (13, ip . With the objective of applying Lemma 13.11 we need to bound the 
number of sets in 7 V = { F D A v : F S J}. 

Assume first that G is A-regular. Let A v = {v, u±, . . . , «a}- We will rule out A + 1 
candidates for F D A v , namely 

0,{ui},...,{u A } (3.3) 

This then shows that \!J V \ < 2 A+1 — A — 1 and thus the bound follows from Lemma [3TTT , 

To see that (|3.3p holds, observe that F £ J contains a B C F such that the induced 
subgraph G[B] is maximal bipartite. To reach a contradiction, assume that there exists a 
v € V with F C)A V Q {ug}. Since B C F, we have Bni„ C {t^}, implying that u does not 
belong to B, and that at most one of its neighbours does. Consequently, G[B U {v}] is also 
bipartite, and v belongs to a partite set opposite to any of its neighbours. This contradicts 
the fact that G[B] is maximal bipartite. 

To establish the non-regular case, we can proceed as in the proof of Theorem 13.21 
adding each u € V with d(u) < A to some A — d(u) of the sets A v not already including 
u. Note that by adding y new vertices to A v originally containing x vertices, we get 
\{F n A v : F e 7}\ < 2y(2 x - x - 1). Next, since 2» {2 X -x-l)< 2^ +:r - (y + x) - 1 for all 
non-negative integers y, x and log {2 X — x — 1) is a concave function, the bound follows as 
before via Jensen's inequality. ■ 

Theorem 3.9. For an n-vertex graph G with maximum degree A we can decide in time 

0((2 A+1 -A-l) n /( A+1 )) 
whether G admits a covering with k independent sets. 

Proof. When k is even, it is easy to see that G can be covered by k independent sets if and 
only if it can be covered by k' = k/2 maximal bipartite sets, so we will use Algorithm C on 
max 93 and investigate whether c(V) ^ 0. 

When k is odd, we again use Algorithm C with k! = {k — l)/2 maximal bipartite sets, 
but this time we check whether an X is output such that both c{X) ^ and V \ X is 
independent in G. 

In both cases the running time bound follows from Theorem 13.81 ■ 

4. Concluding Remarks 

Since the presented improvements on running time bounds are modest, one can ask 
whether this is because of weak bounds or because of inherent limitations of the technique. 
We observe that the running time bounds in Theorems 13.41 [3T6l and !3.9l are met by a disjoint 
union of complete graphs of order A + 1. Thus, either further trimming or splitting into 
connected components is required for improved algorithms in this context. 

We chose to demonstrate the technique for Chromatic and Domatic Number since these 
are well-known and well-studied. To briefly demonstrate some further application potential, 
more artificial problem variants such as determining if a A-regular graph has domatic 
number at least A/2, or if the square of a A-regular graph has chromatic number at most 
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3A/2, admit stronger bounds. For example, if G has domatic number at least d, d even, 
then its vertices can be partitioned into two sets, both of which contain d/2 dominating 
sets. This suggests the following meet-in-the-middle strategy. Run Algorithm D with 3 
equal to all dominating sets and k = d/2, but modify Step D7 to insert X U {j} only 
if \A V \ (X U {j})\ > d/2 holds for all vertices v. At termination, we check whether the 
algorithm has output two sets X and Y such that X UY = V and d(X), d(Y) > 0. (For 
example, one can check for duplicates in a table with entry {X, V\X} for each output X with 
d(X) > 0.) This algorithm variant considers only sets with many forbidden intersections 
with the neighbourhoods of vertices, which translates into stronger bounds via Lemma 13. 11 
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